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Abstract —The main objective of this paper is to propose 
an alternative procedure to carry out one of the key steps of 
immersion and invariance stabilising controller design. Namely, 
the one that ensures attractivity of the manifold whose internal 
dynamics contains a copy of the desired system behaviour. 
Towards this end we invoke contraction theory principles and 
ensure the attractivity of the manifold rendering it horizontally 
contractive. The main advantage of adopting this alternative 
approach is to make more systematic the last step of the 
design with more explicit degrees of freedom to accomplish 
the task. The classical case of systems in feedback form is used 
to illustrate the proposed controller design. 

Index Terms —Stabilization; contraction; nonlinear sys¬ 
tems. 


I. Introduction 

Immersion and invariance (I&I) is a controller design 
technique that has been recently proposed in the literature 
to stabilise non-linear systems [1]—see also the recent book 
[2] where many applications of the technique are presented. 
The I&I approach captures the desired behaviour of the 
system to be controlled by introducing a target dynamical 
system. Then, a suitable stabilizing control law is designed to 
guarantee that the controlled system asymptotically behaves 
like the target system. More precisely, the I&I methodology 
relies on finding a manifold in the plants state-space that 
can be rendered invariant and attractive by feedback control, 
such that (i) on the manifold, the closed loop dynamics be¬ 
haves like the desired dynamics (ii) away from the manifold, 
the control law steers the state of the system towards the 
manifold. The usual way to carry out the latter step is to 
define an extended dynamical system given by a copy of 
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the plant and by a new error dynamics, denoted by the 
coordinate z, that measures the distance to the manifold. 
Then, a full-state feedback controller must be designed to 
ensure boundedness of the plant state and convergence to 
zero of the z coordinate. The main stabilisation result in 
I&I states that the evaluation of this control law on the 
manifold defines an asymptotically stabilising controller for 
the system. The construction leads to a static controller, since 
the control law is a function only of the plant state. 

The design of the aforementioned full-state feedback 
controller is not systematic and finding a controller that 
renders the desired manifold attractive could be challenging 
in practice. The main objective of this paper is to carry 
out this step by exploiting contraction theory principles [6]. 
More precisely, we will use horizontal contraction [4] to 
draw geometric conditions that guarantee the attractiveness 
of the desired manifold. The main advantage of adopting this 
alternative approach in I&I is to make more systematic the 
design of the control action away from the desired manifold. 
We anticipate that the stabilization of the extended system of 
I&I is replaced by the stabilization of the prolonged system 
[3], defined by the plants and its linearization. In comparison 
to I&I, the local nature of the approach pursued in this 
paper provides more degrees of freedom in the design of the 
controller, possibly widening the use of I&I in applications. 

The paper is organized as follows. Section HQ briefly 
recalls the standard I&I controller design procedure. The 
novel design based on horizontal contraction is illustrated 
in Section [III] The section provides the main result of the 
paper, whose proof is in Appendix B. The classical example 
of systems in feedback form is presented in Section [IV] 
Concluding remarks are detailed in Section [V] Appendix A 
contains a counterexample to the classical I&I design. This 
minor issue of the classical I&I is easily fixed by enforcing 
a simple extra assumption. 

Notation For x £ !!" we denote the Euclidean norm \x\ 2 := 
x T x. Given a function / : M" -A R we define the differential 
operators 


V/ := 






where Xi £ W‘ is an element of the vector x. For a mapping 
g : R n —> R m , its Jacobian matrix is defined as 


Vff := 


(V 5 i) T 

(V.g m ) T 




where gi : R ra —>• R is the i-th element of g. When clear 
from the context the subindex of the operator V and the 
arguments of the functions will be omitted. All the functions 
in the paper are assumed sufficiently smooth. 

II. The Standard I&I Stabilisation Procedure 
Consider the system 

x = f{x)+g{x)u (1) 

with state x £ R", the input control u £ R m , and an 
assignable equilibrium point 

£* £ {x £ R n | g ± (x)f(x) = 0} 

to be stabilized, where g 2 - : R™ — > R("- m ) xm i s a full- 
rank left annihilator of g{x). Stabilisation is achieved in I&I 
fulfilling the following four steps. The reader is referred to 
[2] for the proof of the next proposition 

Proposition 1: Assume that there exist mappings 

a:R p H>R p , tt : W -> R" , c : R" —>• R m , 

</> : R n R"-p , v :R n x R”- p R m , 

with p < n, such that the following hold. 

(Al) (Target system) The system 

i = a(0 , (2) 

has a globally asymptotically stable equilibrium at & £ 

R p and x * = 7t(£*). 

(A2) (Manifold invariance condition) For all £ £ R p , 

/W0) + fi r ( 7r (0) c ( 7r (0) = V7r(£)a(£) ■ (3) 

(A3) (Implicit manifold description) The following set iden¬ 
tity holds 

JA := {x £ R"|a; = 7r(£)} = {x £ R"|</>(a;) = 0} . 

(4) 

(A4) (Manifold attractivity and trajectory boundedness) Con¬ 


sider the system 

x = f(x)+g(x)v(x,z) (5) 

z = Vf(x)[f(x) + g(x)v(x, z)\, (6) 

with the initial condition constraint 

z(0) = 0(z(O)), (7) 

and v(x, z ) verifying 

uM£),0)=c(7t(£)), V£ G R p . (8) 

All trajectories of the system are bounded and satisfy 
lim z(t) = 0 . (9) 

t—foo 

Then, £* is a globally asymptotically stable (GAS) equilib¬ 
rium of the closed-loop system 

x = f(x)+g(x)v(x,<j>(x)). (10) 


Following the discussion in the introduction, the accom¬ 
plishment of step (A4) is not systematic and may challenge 


the successful completion of the I&I design. Exploiting [4], 
we propose in the next section to replace (A4) by a novel 
condition based on horizontal contraction [4]. 

Remark 1: In comparison to the results presented in [1], 
[2], we have added the initial condition constraint i[7Ji and 
the requirement <[8}. The first condition ensures that z(t) = 
tj>(x(t)), Vf > 0, while the second one guarantees that the 
.x-system behaves like the ^-system when restricted to the 
manifold JA. These requirements were implicitly assumed 
in previous works. If these conditions are not imposed it is 
possible to show that the claim of Proposition [I] is false. An 
example of this fact is given in Appendix A. _i 

Remark 2: The I&I technique makes contact with the 
literature of invariant manifolds stability [8] and of condi¬ 
tional stability (relative to a set) [7]. Indeed, in the I&I 
technique the action of a state-feedback controller renders 
invariant and stabilizes a suitable submanifold of the system 
state space while enforcing a desired steady-state behavior, 
represented by the target dynamics. j 

III. The I&I Horizontal Contraction Procedure 

The proposition below proposes to replace the step (A4) 
in Proposition|T|by a horizontal contraction based design that 
ensures attractivity of the manifold JA. 

Proposition 2: Given the conditions (A1)-(A3) in 
Proposition [I] assume there exist mappings 

P : R" -a- r("-p) x (™-p) j p — pT > q 
R : R” —> R(«-P)xn ) 

/3 : R n -a- R m , 

p : W 1 —> R, p > 0 , 

such that the following holds. 

(A4’) (Manifold attractivity via horizontal contraction) 

(i) For all x £ R n , R(x) is full rank and 

R«0) = v^(7r(0), £ IT. (11) 

(ii) For all £ £ W 

/3(tt(0) = c(tt(0). 

(iii) The candidate Finsler-Lyapunov function V : M" x 
R” —y R> 0 given by 

V(x,6x) := Sx T R T (x)P(x)R(x)Sx, (12) 

satisfies 

V(x, Sx) < — p(x)V(x, Sx) (13) 

along the trajectories of the prolonged system 

x = f{x)+g{x)l3{x) (14) 

5x = y[f(x) + g(x)/3(x)]8x. (15) 

(iv) The trajectories of (IT4l > are bounded. 

Then, £* is a GAS equilibrium point of the closed-loop 
system m. Furthermore, if the fixed point & of the target 
system Q is hyperbolic, then re* is hyperbolic. _i 

Remark 3: A natural simple choice for R(x) is V0(:r), 
provided that V<fi(x) is full rank for all x £ R”. _i 


Remark 4: In contrast with the classical I&I Proposi¬ 
tion [2] directly provides the static state-feedback controller 
0(x). This should be compared with the control v(x, z ) that 
should verify condition (A4) for the augmented system (O, 
©, which is later evaluated on the manifold to generate the 
actual control to be applied, that is, v(x, j 

Remark 5: Proposition |3 can be formulated in a similar 
way for any forward invariant region C C R ra . If C = R" 
then, as stated in the proposition, one gets GAS. Otherwise, 
one gets regional stability. This formulation may be useful 
in applications when global results are difficult to achieve or 
when the system lives in a manifold different from R”. Note 
that if C is compact, then the condition (iv) of boundedness 
of trajectories is automatically satisfied. j 

IV. Application to Systems in Feedback Form 


Fixing p(x) = k > 0 the condition (fl3l) is satisfied if and 
only if 

M{x 1 ) + M(xi)[Q(x) + ^-1} + [Q T (x) + ^-I]M(x 1 ) < 0. 

( 22 ) 

We are in position to state the following proposition. 

Proposition 3: Consider a system described by equa¬ 
tions of the form (IT6l > and suppose there exist mappings 
7 t 2 : R" —> R and p : R( n+1 ) —> R such that the following 
holds. 

(a) The system 

±1 = f(xi, 7t 2 (xi)) 

has a GAS equilibrium at zero. 

(b) The inequality d23 is satisfied for some k > 0. 

(c) For all £ £ R 


Consider the class of systems in feedback form de¬ 
scribed by the equations 

ii = f(x 1 , 3 : 2 ), , 16 , 

X 2 = U, 

with x := col(tci, X2) £ R™ x R, and u £ R. Consistent with 
the standard backstepping scenario [5] assume there exists a 
mapping 7r 2 : R” —> R such that the system 

Xl = f{ Xi,TT 2 (Xi)) 


has a GAS equilibrium at the origin. A sensible choice of 
the target dynamics is then given by 

£ = mM0), (i7) 


and this implies that the mapping 7 t(£) has the form 




£ 

MO 


To verify Assumptions (A2) and (A3) of Proposition Q] we 
can choose 


<0M0) = ^MOfiOMO) ( 18 ) 

4>(x) = X2-7T 2 (lEi), (19) 


which clearly satisfy ((3 and (Q]l 

The differential relation of the system (IT6t in closed- 
loop with the control /3(x) is 


5x = 


V ai /(x) V X3 f{x) 
V Xl /3(x) V X2 /3ix ) 


5x =: Q{x)5x , 


From Proposition [2] we select Ii(x) = Vo(x) and P(x) = I. 
Whence the Finsler-Lyapunov function ( 1 1 2b takes the form 


V{x, 5x) = 5 x t M(xi)5x , (20) 


where we have defined 


M(x 1 ) 


V0(j:)] t V0( x) 

V7r 2 (a:i)[V7r2(xi)] T 

-[V7r 2 (xi)] T 


-V7T 2 (xi) 

1 


( 21 ) 


m,MO) = ^M0f(0M0)- 

Then, the system (fTTil i in closed-loop with (3{x) has a GAS 
equilibrium at zero. _i 

Example 1: To illustrate the result in Proposition [3 
consider the two-dimensional system 

xi = -xi + \x\x2 ( 23 ) 

±2 = u, 

in which A > 0. We proceed now to verify condition (a). 

Selecting 7 t 2 (xi) = — x\ we obtain 

xi = f{x i,7T 2 (xi)) = -xi - Ax® 

which has a GAS equilibrium at zero. To check condition 

(b) we, first, compute 

x\ + X 2 
4xf 2xi 

2xi 1 

—1 + 3Axfx2 Axf 

V X1 J3 ix) V x Jix) _ ■ 

Some lengthy, but straightforward calculations, show that 

Pix) = “(A: - 4)xf - ifcx 2 - 2Axfx 2 (24) 

solves (l22l) with identity. It only remains to verify condition 

(c) , which holds true because 

PiOMO) = 2C 2 (1 + AC 4 ) = (-2Q J-( + A( 3 (-( 2 )] . 

V7r 2 (C) f(£,n 2 ( 0 ) 

In conclusion, the system (123 in closed-loop with the 
control (123 has a GAS equilibrium at the origin. _i 

V. Conclusions 

An alternative procedure to complete the design of I&I 
controllers for stabilization of nonlinear systems has been 
proposed. The central idea is to replace by a contraction- 
based design the stabilization step on the extended dynamics 
©,© required by condition (A4) of the I&I procedure. 
The main advantage of the contraction-based approach is to 
render more systematic the design and to give more degrees 


(fix) = 
Mix) = 

Qix) = 












of freedom for its accomplishment. The key step of the novel 
design is the use of horizontal Finsler-Lyapunov functions 
[4] that decays along the trajectories of the prolonged system, 
in the spirit of classical Lyapunov theory. Of course, similarly 
to all constructive procedures for the design of nonlinear 
controllers or observers, for the successful application of the 
novel design proposed by the paper it is necessary to solve 
a partial differential equation. In particular, for systems in 
feedback form, it is necessary to find a controller /3(x) that 
satisfies d22l> (f3(x) is encoded in Q{x)) for a suitable choice 
of R{x) and P. 

From the conceptual viewpoint, the use of Finsler- 
Lyapunov functions replaces the stabilization of the ol'f- 
manifold coordinate z of I&I with the horizontal stabilization 
of the linearization along trajectories. For instance, the 
method proposed in this paper stabilizes the linearization of 
the system along suitable directions of its tangent space, thus 
providing a local and intrinsic feedback design procedure 
that does not require any a-priori definition of the ol'f- 
manifold coordinate z. The advantage is a more general 
design method, possibly. This generality is directly encoded 
into the conditions of Proposition [2j the z coordinate of clas¬ 
sical I&I is replaced at local level by the matrix R(x), which 
is one of the free parameters to be selected in the formulation 
of the partial differential equation (l22l) (Al (x) depends on 
R(x)). The intrinsic nature of the design combined with the 
increased degrees of freedom make the present formulation 
of horizontal contraction-based I&I a promising stabilization 
tool for applications. 


where 9 f 0 is a constant parameter and u = col(ui, U 2 ) 
is the control input. The control objective is to stabilize the 
system at the origin using the I&I procedure. 

First, we select the target system as f = — £, which 
clearly has a GAS equilibrium at zero, verifying the first 
part of (Al). Selecting 7 t(£) = col(£,0) it is easy to see 
that the manifold invariance condition (A2) holds with the 
constant control 

c(tt(0) = Q ■ 

Moreover, 7r(0) = col(0,0), verifying the second part of 
(Al). The implicit manifold condition (A3) is verified with 
t f>(x ) = X 2 ■ Finally, we need to define a controller v(x,z) 
for the augmented system 

±1 

X2 
z 

that ensures boundedness of trajectories and lim^oo z(t) = 
0. This is clearly guaranteed with the selection 

/ x r ° 1 

V{x,z)= . 


—xi + 9 + vi(x, z ) 

V 2 (x,z) 

—x\ + 9 + vi(x, z ) 
V2 (x,z) 


= V<p(x) 


= v 2 (x,z), 


It is claimed in Theorem 2.1 of [2] that applying the 
control 


U = v(x, (j>{x)) 


0 

-X2 
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Appendix 


to the x-system ensures the origin is a GAS equilibrium. 
But the resulting closed-loop system 

x\ = —x\ + 9 

X2 = —X2 

has a GAS equilibrium at (0,0), not at the origin. 

The source of the problem is that, if we do not impose in 
(A4) the initial condition (IT} we have only that z = <j> = x 2 
but z(t) ^ <j>(x(t)) = X2(t). Indeed, integrating the system 

x2 = —z 

z = —z, 

we get 

z{t) = e-*z( 0) 

X 2 (t) = 2 2 ( 0 ) - f* e~ T z( 0 )dr = x 2 ( 0 ) - z(0)(l - e _t ). 

(25) 

Notice that the condition © is also not verified since 


A. Counterexample to Theorem 2.1 of [2] 

Our objective in this appendix is to show that if we 
follow all the steps of the standard I&I procedure of Propo¬ 
sition G ]— without imposing the conditions 0 and ([8}— we 
cannot guarantee GAS of the equilibrium. Towards this end, 
consider the two-dimensional, linear, time-invariant system 

21 = —2i + 9 + U 1 

22 = U 2 


2(tt(£),0) 


^ c(tt(£)) 


-9 

0 


It is clear from (l25l > that imposing the initial condition 
©. that is, z( 0 ) = 22(0) we get z(t) = 22(f). But we still 
need to modify the controller to comply with ®. A simple 
choice being 

, s [-01 

v(x,z)= . 
















For this new choice the control 

u = v{x , (f>(x)) = 

yields x = —x, which certainly has a GAS equilibrium at 
zero. 



B. Proof of Proposition \2\ 

The proof is divided in four parts establishing global 
attractivity of (I) the manifold and (II) the equilibrium point, 
(III) local stability and (IV) hyperbolicity of the equilibrium 
point. 


I. Global attractiveness of M 

Take \8x\ x := y/V(x, Sx). Given any differentiable 
curve 7 : [0, 1 ] —> R n define the horizontal length £( 7 ) := 
fo Note that £( 7 ) 0 iff R('y(s))j(s) 0 0 for 

some s £ [0,1]. Thus, £( 7 ) 0 if 7 (s) A4 for some 

s £ [0,1]. For instance, consider any x and y in At. By 
construction, there exists a differentiable curve 7 such that 
7 ( 0 ) = x, 7 ( 1 ) = y and 7 (s) £ M for all s £ [0,1]. Then 
f?( 7 (s)) 7 (s) = V^( 7 (s)) 7 (s) = 0 , by (Qj} , which implies 
It(' s ) l-y(s) = 0, thus £( 7 ) = 0. In a similar way, consider any 
x £ M. and y f A4 and let 7 be any differentiable curve 
such that 7 ( 0 ) = x and 7 ( 1 ) = y. Then, because of the rank 
condition on R(x) and differentiability of //(:;;), there exists 
a measurable subset of I C [0,1] such that | 7 (s)| 7 ( s ) 0 0. 
Thus, £( 7 ) > 0. 

Let ipt(xo) denotes the flow of the system x = f(x) + 
g(x)/3(x) at time t from the initial condition ipo(x 0 ) = 
.To £ K n . Exploiting the boundedness of trajectories, 
global attractiveness of M can be proven by showing that 
lim^oo £(" 0 t ( 7 )) = 0 for any given curve 7 such that 
£( 7 ) 0 0. We show this in the next two paragraphs. 

By boundedness of trajectories, for any 7 : [0,1] -> R" 
there exists a compact set JC such that, ^ 4 ( 7 ( 5 )) £ K. for 
each s £ [0,1] and t, > 0. By continuity, for each x £ K. and 
sx £ r". m guarantees that there exists A > 0 such that 
V(x, Sx) < —XV(x, Sx) < 0. It follows that 


ds 


V tMTOO), izMli 8 )) < exp(-A£)V j(s), 7-7(5) 


ds 


which implies that 




. A t 

< exp —— 


V’tGCs)) 




7(x) 


Thus, £(ip t ( 7)) < exp(-^)£( 7 ). 

Suppose now that 7(0) £ A 4 and 7(1) A 4 . By (A 2 ), 

t/>t(7(0)) £ M for all t > 0 (manifold invariance). Thus, 
the combination of lim^oo £(1/7(7)) = 0 with boundedness 
of trajectories guarantees that "0* (7(1)) converges asymptot¬ 
ically to M. 


II. Global attractiveness of x * 

By (A2), (A4’) and the boundedness of trajectories we 
have that any trajectory of the closed loop system converges 
to the manifold (j>{x) = 0. Moreover, by (Al) and (A2), 
the manifold is invariant and internally asymptotically stable, 
hence all trajectories of the closed loop system converge to 
the equilibrium 2*0 

III. Local stability of x* 

To conclude the proof we need to show that x * is 
Lyapunov stable. Note that any trajectory of the closed loop 
system is the image through the mapping 7r(-) of a trajectory 
of the target system £. Moreover, for any e \, there exists 
S\ such that |£(0)| < <5i implies |£(f)| < e±. Hence, by 
regularity of n, for any e > 0 there exists 6 > 0 such that 
K(£( 0 ))[ < S implies |tt(^(£))| < e. 

IV. Hyperbolicity of 

Without loss of generality take a;* = 0 and define S := 
Va(0), n := Vtt( 0), and A := V[/(0) + s(0)c(0)]. For 
semplicity denote the matrix R( 0) by R and P(0) by P. 
The linearization of the closed loop system computed on the 
fixed point reads x = Ax. £ = S'£ denotes the linearization 
of the target system at the fixed point £* = 0. 

Now, S is Hurwitz by assumption. The span of the 
columns of n and R 1 define two orthogonal subspaces of 
the state space. To see this, note that 11£ £ T Xr M therefore 
iffl£ = V</)(0)n^ = 0. Since £ has dimension p, it follows 
that Rl I = 0. Clearly, n T P T = (fffl) T = 0. It follows 
that the state of the linearized closed loop system can be 
decomposed as 

x = n£ + R r e 

where e is a vector in E f "~C. In particular, take 

e := {RR T )~ 1 R(x - n|) = {RR T )~ 1 Rx . 

The invertibility of RR T follows from the rank condition on 
R. The last identity follows from orthogonality, i.e. f?n = 0. 

Take any trajectory xf). Then, there exist positive 
constants ci, C 2 , C 3 , C 4 , C 5 , cq such that 


x{t)\ < 

ci|n£(f) + R T e(t) | 

< 

Ci|n£(f)| + Ci|P T e(f)| 

< 

c 1 \ui(t)\+c 1 \(RR T )~ 1 Rx{t)\ 

< 

c il n ?WI +c 2 y/V{ 0 ,x(t)) 

< 

C 3 |C(*)I +c 4 exp(-^)v / V( 0 ,i( 0 )) 

< 

c 3 exp(A max (S)t) £( 0 )| +c 5 exp(-^ 

< 

c 6 exp(max{A max (S'), -^}£)|i( 0 )| 


where A ma x(< 5 ') is the largest eigenvalue of S and A is the 
local decay rate of V (part 1 of the proof). Finally, expo¬ 
nential stability of the linearization implies local exponential 
stability of x * for the closed-loop dynamics. 


*This steps coincide with the ones of Theorem 2.1 of [2], The same 
remark applies to the derivations of part III. 











